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SYSTEM OF PHASE OSCILLATORS WITH DIAGONALIZABLE 

INTERACTION 

TAKASHI NISHIKAWAt AND FRANK C. HOPPENSTEADTt 

Abstract. We consider a system of N phase oscillators having randomly distributed natural 
frequencies and diagonalizable interactions among the oscillators. We show that, in the limit of 
N — » oo, all solutions of such a system are incoherent with probability one for any strength of 
coupling, which implies that there is no sharp transition from incoherence to coherence as the coupling 
strength is increased, in striking contrast to Kuramoto's (special) oscillator system. 
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1. Introduction. Synchronization of coupled oscillators is a ubiquitous phe- 
nomenon in natural and artificial systems. Examples include synchronization of pace- 
maker cells of the heart 112) . rhythmic activities in the brain 0JEj) synchronous 
flashing of fireflies PE], arrays of lasers [El El, and superconducting Josephson junc- 
tions [181 119) . Characterization of the phenomenon using mathematical models has 
been a topic of great interest for researchers in various scientific and engineering dis- 
ciplines. 

Wiener 117) . who recognized the ubiquity of synchronization phenomena in 
the real world, made a first attempt at characterization using the Fourier integrals. 
A more successful approach was taken by Winfree [20], who used a population of 
interacting limit-cycle oscillators to describe synchronization properties. He realized 
that if the interactions among the oscillators are weak and the oscillators are nearly 
identical, the separation of fast and slow timescales leads to a reduced model that 
can be expressed in terms solely of the phase of each oscillator. Kuramoto ^Tj] put 
this idea on a firmer foundation by employing a perturbation method to show that 
the reduced equation has a universal form. His analysis of this model in the case of 
mean-field coupling kicked off an avalanche of theoretical investigations of his model 
and its generalizations. 

More generally and rigorously, if each oscillator has an exponentially stable limit- 
cycle and interactions among them are weak, the reduced phase equation can be shown 
(see Theorem 9.1 H3 P- 253]) to have the form 

6 = oj + ef(0), 8eT N , (1.1) 

where oj € is the vector of natural frequencies of the oscillators that are coupled 
to one another through the interaction function / : T N — > M. N , and e > represents 
the overall strength of the coupling. The universal form of the interaction function, 
derived by Kuramoto |10| under the additional assumption that the oscillators are 
almost identical, corresponds to the choice 

JV 

/i(0) = X)M*J-*O. i = l,...,N, (1.2) 
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where f(9) — (fi(0), . . . , /tv($)) T - The mean- field model that he studied results when 
hij(x) = sin(x)/N for all 

Let 9(t) be a solution of (|i.l() . The oscillators i and j are said to be locked if 
lim^oo 0i(t)/9j(t) = 1. The solution is said to be coherent if all pairs of oscillators 
are locked. If none of the oscillator pairs are locked, the solution is incoherent. A 
solution that is neither coherent nor incoherent is called partially coherent. The main 
conclusion of Kuramoto's work 10 on his mean-field model is that in the limit of 
N — > oo there exists a critical coupling strength e c such that for e < e c the solution is 
incoherent, but for e > s c partially coherent solutions appear, for which the fraction 
of locked oscillator pairs is nonzero. Although his result was important, since this be- 
havior closely resembles the phase transition phenomena widely observed in statistical 
physics, his analysis is heuristic and makes assumptions about the symmetry of the 
distribution of natural frequencies, which might not be necessary for the results |14| . 

In this paper, we consider a class of diagonalizable interaction functions, in which 
separation of variables is possible after an appropriate coordinate transformation. 
This allows us to prove rigorously that for the system 1)1. If) with a generic diago- 
nalizable interaction function, if the solution is partially coherent, then it is almost 
surely coherent. This, together with the fact that the probability of having a coherent 
solution goes to zero in the limit of N — > oo, leads to our main conclusion. Namely, 
for any e > 0, the solution is almost surely incoherent in the limit of N — ► oo. Our 
result shows that a diagonalizable system of phase oscillators cannot exhibit a sudden 
transition from incoherence to coherence, in sharp contrast to the mean-field model 
of Kuramoto. This implies that for the system ill. If) to exhibit a phase transition, the 
interaction function / cannot be diagonalized. 

There is an alternative rigorous approach to Kuramoto's mean-field model, in 
which the partial differential equation for the density of oscillators with certain fre- 
quency, which is obtained by taking the continuum limit N — > oo, is studied to analyze 
the stability of the solutions. See ^H] for an excellent review in this direction. 

The approach taken here is similar to that in |SJ p. 80]. However, some conclusions 
made there might be misleading or lack detailed analysis. This paper is intended to 
correct and clarify those points. 

The rest of the paper is organized as following. In SJ21 we introduce an appropriate 
change of variables to separate a time-like variable from the rest of the system. In Sj3J 
we define diagonalizable interaction and show how complete separation of variables 
can be achieved. We also establish some properties of diagonalizable systems. Then, 
in we introduce randomness of the natural frequencies of the oscillators and state 
our main results. Finally, we discuss some approximate behavior of the system for 
large N in ^ and SQis reserved for concluding remarks. 

2. Separation of the Time-like Variable. In this and the following sections, 
we consider the system l|l.l[) of N phase oscillators, where the natural frequency vector 
u> and the coupling strength e are fixed (nonrandom) constants. We will consider oj to 
be a random vector in 21 in order to make probabilistic statements about the system. 
Let us suppose that the interaction function / satisfies two conditions, 

(CI) l T f(6) = for all 9 e T N and 

(C2) f(vl + 9) = f{9) for all 9 € T w , 
where 1 = (1/^/N, . . . , l/\fN) T . The condition (CI) says that the interaction func- 
tion is orthogonal to the vector 1. The second condition (C2) expresses the translation 
invariance of / along the direction of 1. If, for example, the interaction function has 
the form l|1.2[> . these conditions are satisfied if the functions hij are odd. In particular, 
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the mean- field model of Kuramoto does satisfy these conditions. 

Under conditions (CI) and (C2), the system Ijl.lfl can be separated into two inde- 
pendent systems — one for the time- like variable and the other for the phase deviations. 

Let W be an N x (N— 1) matrix whose columns, denoted by Wj , j = 1, . . . , N — 1, 
form an orthonormal basis of the subspace l -1 = {x S K : l T x — 0}. In other words, 
W is an N x (N — 1) matrix that satisfies 1 T W = and W T W — /jv-i, where Jjv-i 
is the (N — 1) x (N — 1) identity matrix. Then, the change of variable 

9 = vl + Wu (2.1) 

converts the system ljl.ll) into two systems, 

v = l T w, (2.2) 
u = W T uj + eW T f(Wu), (2.3) 

which can be solved separately. 

Systems satisfying the conditions (CI) and (C2) arise in mathematical neuro- 
science El , in which 9 often takes the form uit + tfi in the limit t — > oo, where 
ui is the vector of carrier frequencies and </> is the vector of phase deviations. The 
equation (|2.3|) . in some sense, governs the behavior of the phase deviations. 

The solution to (|2.2|l is v(t) — v(0) + (l T ui)t, and hence the variable v is time-like 
if 1 t uj or, equivalently, if the average natural frequency ^i^i/N is nonzero. 
Thus, the behavior of the solution of is essentially determined by l|2.3|l . 

Recall that the solution is called coherent if lim^oo 6i(t)/6j(t) — 1. This can 
be rephrased in terms of the vector fi = limt^oo u{t)/t of output frequencies of the 
u-equation l|2.3|l . if it exists. 

Lemma 2.1. Let u(t) be a solution of (|2.3(l . and suppose that = limt^oo u(t)/t 
exists. Then, the solution of (|1 . 1|> is coherent if and only if fi = 0. 

Proof. Let fl = Iim t -oo 0(f) A = (1 T ^)1 + WfJ,. Then fj, = implies that 
O = (l T ui)l, which in turn implies that lim^oo 9i(t)/9j(t) = 1. 

Conversely, if limj—,^ 6i(t)/9j(t) = 1, Cl must be a multiple of 1. Since W is 
orthogonal to 1, this implies that Wfi = 0. Since W is invertible, it follows that 
fi = 0. □ 

As an immediate consequence of Lemma 12.11 if the solution of (|2.3(l tends to an 
equilibrium, then the corresponding solution of (|l.l(l is coherent. For example, if 
the interaction function / that satisfies (CI) and (C2) is a gradient vector field, i.e., 
f(9) — — Wq(9) for some potential function Vq : T N — > R, then the w-equation l|2.3(l 
is also a gradient system: 

u = -V [~lu t Wu + eV{u)] , (2.4) 

where V(u) = Vq(Wu). A minimum u* of the potential function — ui T Wu+eV(u) then 
corresponds to the vector of phase deviations for a coherent solution of the original 
oscillator system Q1.1JI . As in the proof of Lemma \2. II we have $1 = (l T w)l in this 
case, meaning that the output frequency of every oscillator tends to the mean natural 
frequency uj = J^- u>i/N of the oscillators. For the interaction of the form (|1.2fl . the 
potential function takes the form 

1 JV JV N N /N-1 

V ^ = E E - ^) = -^v EE H v E - w J»] u >< 

i=l j = l i=l j=l \k=l 

where H. tj (x) = J* h lj (y)dy. 
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3. Diagonalizable interaction. Assuming that the interactions among the os- 
cillators are diagonalizable enables us to carry out a rigorous analysis of the system. 

Definition 3.1. We say that the system (or the interaction function f) is 

diagonalizable if there exist an N x (N — 1) matrix W and real, continuous, periodic 
functions pj such that 

(i) 1 T W = 0, 

(ii) W T W = Jjv-i, and 

(iii) f(Wu) = Wp{u) with p(u) = (pi(ui),. . . ,p N ~i(u N -i)) T . 
For example, W = defined by 



= -= ( sin + cos — ^- I 



1 ( . 

—= sin 

Vn \ 


2irjk 


N + 


2 • ( 
—= sin 

Vn V 


2irjk 


N + 



2irjk\ 



(3.1) 



satisfies these conditions. 

When the system (|1 . 1|> is diagonalizable, the equations for the components of u 
become independent of other components: 

Uj = aj+epj(uj), j = 1, . . . ,N - 1, (3.2) 

where we set aj = WJlo. Thus, the problem is reduced to solving a scalar differential 
equation for each j. The following lemma applies to each equation in (|3.2() . 

Lemma 3.2. Let e > 0, and let a be a real number. Letp(u) be a real, continuous, 

periodic function with period L > 0. Define m = min p(u), M = max p(u). For 

o<u<l y ' o<u<l y ' 

any solution u(t) of u = a + ep(u), the limit fj, p (a,s) = lim^oo u(t)/t exists and 

i \ )L/T(a,e), a < -eM, a > -em, 
fi p (a,e) = i (3.3) 
0, — eM < a < —em, 



whe 



T(a,e) 



du 



o a + ep(u) 



is the "period" of the solution in the case of a < —eM or a > —em, in the sense that 
u(t + T(a,e)) =u(t)+L. 

Proof. If —eM < a < —em, then any solution u(t) tends to a zero of the function 
a + ep{u). Hence, fip(a, e) = 0. 

For notational simplicity, let us drop the dependence of T(a, e) on a and e below. 
Suppose a > —em, so that a + ep(u) > for all u. It is straightforward to show that 
the function u(t) defined implicitly by the formula 

du 

= t 



a + ep{u) 



is the unique solution of u = a + ep(u) with the initial condition u(0) = uq, and that 
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it satisfies u(t + T) = u(t) + L. We have 



fi p (a,e) = lim u(t)/t 

t—>oo 



lim — un 

t^oo t 1 



a + e lim — 

t— too t 



[a + ep{u{s))]ds 



p{u{s))ds. 



Let n be the largest integer for which nT < t. Then, by changing the variables in 
each integral using the translation by multiples of T, we see that 



/ p(u(s))ds = V / p 

JO k=l J(k~l)T 



(u(s))ds - 



p(u(s))ds 



nT 



n / p(u(s))ds - 



t-nT 



p(u(s))ds. 



Consequently, 



p(u(s))ds — 



T 



p(u(s))ds 



< 



( n 


1 


u 


~ T 


nT 


-t\ 



p{u{s))ds + - 



t-nT 



p{u{s))ds 



tT 



\p(u(s))\ds 



t-nT 



\p{u{s))\ds 



< Tmax{|m|, \M\} + I Tmax{|m|, \M\} 



0. 



as t 
i r T 



o, showing that the limit lim^oo j J p(u(s))ds exists and is equal to 
Jo p(u(s))ds. Thus, by changing variables from s to u and translating by Uq, 



we see that fi p (a, e) exists and 

u p (a,e) 



e 

TJu 

£ 



p{u{s))ds 

p{u)du 
T J a + ep(u) 



= L/T. 

If a < —eM, then, by replacing u with —u, a with —a, e with — e, and m with 
M, the problem reduces to the previous case. The lemma is proved. □ 

The function e) in Lemma 13.21 which can easily be shown to be differentiable 
with positive derivative outside the interval [— eM, eM], determines the relationship 
between the input frequency a and the output frequency u p (a,e). If we take p(u) = 
sin(u), for example, the integration in the expression of n p can be carried out, and 
we get 



-%/a 2 — e 2 , a < —e, 



fi P (a,e) = < 0, 



~e < a < e, 



V a 2 — e 2 , a > e, 
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Fig. 3.1. (a) The graph of the input-output frequency function /i p (a,e) vs a for s = 1 and 
p(u) = sin(w). (b) The corresponding density g(fj,;l) when a is the standard Gaussian random 
variable. 



the graph of which is given in Fig. 13. If a) for e = 1. With this function, the output 
frequency vector \i of the u-equation l|2.3[) can be written as 

H = n(a) = (/i Pl (ai,e), . . . , n PN _ 1 (a N ^ 1 , e)) T . 

It is important to note here that \i does not depend on the initial condition 
u(0) = uq, which implies that it is also independent of the initial condition for 9. In 
other words, the initial condition for the system (f 1 - X|) does not affect the behavior of 
its solution, as far as its coherence properties are concerned. Therefore, in this sense, 
coherence, partial coherence and incoherence are properties of the system rather than 
of individual solutions for a diagonalizable system. 

4. Randomly Distributed Frequencies. In this section we consider w to be a 
random vector in M. N . We take the components u>x, . . . , gjn of u to be independent and 
identically distributed (i.i.d.) random variables with mean and variance a 2 > 0. In 
the general case of mean uiq ^ 0, the problem can always be reduced to the zero-mean 
case by the translation of 9 by — toot. 

Since ui is random, the vectors a and \i are also random vectors in IR^ -1 . Lemma 
13.21 along with the relation a = W t uj can be used to determine the distribution of 
fj, from the distribution of lu. For example, if each uij is standard Gaussian, then so 
is each a,-, in which case the density g(ii;e) for the random variable /J>p(aj, 1) when 
p(u) = sin(u) can be computed. The result is 

, . |(i|e-<»' +s '' /3 „ , J c \ 

9( " ie)= ^pr?? +Werf b)' 

where 8(fi) is Dirac's delta function. The graph of this density is shown in Fig. 13.11 
for e = 1. 

Our main goal is this section is to compute the probabilities that the system (ll.lfl 
is coherent, partially coherent, or incoherent. The following theorem reveals a curious 
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property of a generic diagonalizable system of phase oscillators. 

Theorem 4.1. Let the natural frequency vector uj be a random vector in M. N , 
whose components are i.i.d. with a common continuous distribution. Suppose that i|l.l[l 
is a diagonalizable system of N phase oscillators such that W satisfies the condition 
that Wki ^ Wkj for all k — 1,2, ... ,N and for all i,j = 1,2, ... ,N — 1 such that 
i =/= j . Then, the partial coherence of the system almost surely implies coherence; i.e., 
given that the system is partially coherent, the probability that it is coherent is one. 

Proof. Once again, let Q(ui) — lim^oo 0(t)/t — (l T w)l + W[i(ui). Let S c be the 
set of ui in M. N that corresponds to coherent systems, i.e., S c = {ui £ R N : Q(u) = 0}. 
By Lemma ETT1 we may also write S c = {ui e R N : (j,(ui) = 0}. Let S pc be the set 
corresponding to partially coherent systems, that is, S pc — {u £ R N : £li(ui) — Qj(u>) 
for some i ^ j}. It is easy to see that we can also rewrite this in terms of [i as 



The probability that the system is coherent, given that the system is partially 
coherent, is P(S C )/ P(S pc ) since S c C S pc . This probability is one if and only if 
P(S pc \ S c ) = 0, which would be satisfied if P(S P c \ S c ) — for every pair i ^ j. 
We shall show this next. 

Let us fix i and j. For any A C {1, 2, . . . , N — 1}, denote by Zk the subspace 
{fi e : [ik — 0}, and let Za = U k eA and Z' A — f] k ^ A Zk- Let R A denote the 

subspace {/i e : EkeA( W ki ~ W kj )pi k = 0}. Define Q A = Ra n Z' A \ {0}. We 

will show that P(Qa) = for any choice of A. P(Spc \ S c ) = follows from this by 
taking A = {1, 2, . . . , N - 1}. 

We shall prove P(Qa) = by induction on n — \A\, the cardinality of A. Suppose 
first that n — 1 and, say, A = {1}. Since Wu - W\j ^ 0, we have Ra — S K^ -1 : 

[ii = 0} = Z x and Z' A = f]^ Z k . Thus, Q A = f)k=i z k \ {0} = 0, which implies 
P(Qa) = 0. The same holds for any other A with \A\ = 1. 

Suppose that P(Qa) = for any A with \A\ = n — 1, and consider the case 
\A\ = n. We have 



We see that P(Qa H Zk) = for each k £ A, by the induction hypothesis, since we 
can write Q A n Z k = R A n Z' A n Z k \ {0} = i? Afc n Z^ fc \ {0} = Q Ak with A fe = A\ {k}, 
for which we have \A k \ = n — 1. Thus, if we can show P(Qa \ Za) = 0, then we are 
done. 

We show P(Qa \ Za) = in three steps. First, since Z' A is an n-dimensional 
subspace and Qa C -Ra H ^ is an (n — l)-dimensional subspace, the n-dimensional 
Lebesgue measure of Qa in Z' A must be zero. 




P{Q A ) - P{Qa n Z A ) + P(Qa \ Z A ) 




< 



P(Qa n Z k ) + P(Q A \ Z A ). 
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Next, note that the conditional probability distribution of /z, given \x £ Z' A , is 
continuous with respect to the Lebesgue measure outside the set Z A . This can be 
seen by noting the following: (1) each component can be written as — fi p .(aj) 
by Lemma 13.21 (2) exists and is differentiable except at the origin, again by 
Lemma l3.2l and (3) the conditional distribution of aj = Wkj^k, given that (J,(ui) € 
Z' A (which is equivalent to —eM^ < J2i < —sn^k for all k £ A), is continuous 

everywhere. 

Finally, combining these two observations, we see that P(Q A \ Z A \ Z' A ) = 0, 
which implies that P{Q A \ Z A ) = P(Z' A )P(Q A \ Z A \ Z' A ) = 0. This completes the 
proof of the theorem. □ 

We next describe the behavior of a generic diagonalizable system in the limit of 
N — ► 00. In order to formalize the process of taking the limit, we need to choose a 
sequence of systems of the form (f 1 . 1 |l . Such a sequence can be characterized by the 
following. 

1. Consider a sequence {W^ n '}n=i,2,... of matrices with the following properties: 

a Each WW is an N x (JV — 1) matrix with orthonormal columns. 

(b) l^W = for all N. 

(c) Each W( N ) satisfies the condition for W in Theorem 14. II 

(d) ||VF (JV) ||oo -> as N -► 00. (Here || • IU denotes the maximum matrix 
norm defined by ||^4||oo = max|Aj_y|, where the maximum is taken over 
all elements of A.) This is like a mixing condition that will be necessary 
later in order to apply Proposition 14. 31 

2. Consider a sequence {pj}j=i,2,.., of real, continuous, periodic functions such 

that the corresponding sequence of norms \ \pj\ \ = max\pj(u)\ is bounded. 

3 

3. Consider a sequence {wj}j=i,2,... of i.i.d. random variables with mean loq and 
variance a 2 . 

The sequence of matrices defined by l|3.1|l satisfies the conditions above. Given 

such sequences, for each e > and N, we define Sn,e to be the diagonalizable sys- 
tem (11.1(1 of phase oscillators using the natural frequency vector oj — (oji, . . . , ojn) T , 
the functions {p±, . . . ,pjv-i} and the matrix W^ N K We are now ready to state and 
prove our main theorem. 

Theorem 4.2. Let Sn,c be defined as above. Then, for any fixed e > 0, Sn iE is 
almost surely incoherent as N — > 00; i.e., the probability that Sn,c * s incoherent tends 
to one in the limit of N — > 00 . 

Proof. As mentioned before, we may assume loq = without loss of generality, 
since the ujq ^ case can always be reduced to the loq — case. 

From Theorem 14.11 we know that the probability that 5/v i£ is not incoherent is 
equal to the probability q c that it is coherent. We need to show that q c — ■> as 
N — ► 00. 

Let No < N be fixed. From Lemma f2.ll it follows that q c = P(n = 0). Since 
the sequence {||pj||}j=i,2,... is bounded, we can define M = supMj and m = inf m,j, 
where m* = min Pj(u), Mj — max Pj(u) for each j. Then, Lemma 13.21 implies 

0<u<L 0<u<L 
q c = P(fl = 0) 

= P(-eMj < af ] < -errij, j = 1, . . . ,N - 1) 
< P(-eM < af ] < -em, j = 1, . . . , JV ), 



where a^' 1 = (W- N ^ 
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The following Proposition shows that for each j, converges to a Gaussian 
random variable with mean and variance a 2 . 

Proposition 4.3. Let Xi,X2, ... be a sequence of i.i.d. random variables with 
EXj — and Var(Xj) = E(X?) = a 2 . Suppose that, for each N, real numbers 

&/v,i, ■ • ■ >°n,n satisfy YljLi j = 1- Also, suppose that 



lim max |6jv ?l = 0. 

N->ool<j<N ' J 



Then we have 



N 

S N = J2 b N, 3 X 1 ^M(0,a 2 ) 

3 = 1 

as N — ► oo . 

Proof. Let Y/v.j = bpf,jXj. We will apply the Lindeberg-Feller central limit 
theorem (see, for example, [31 p. 98]) to Y/vj - For this we need to check three condi- 
tions. The first is EY^j — b^jEXj = 0. The second condition is satisfied because 
J2f=i j — Sjli kjv jEX 2 = a 2 > 0. To show that the third condition is satisfied, 
let e > be fixed. We have 



N N 



YnJ>s) =Y,b 2 N ^E(\X J 



where E(X\A) denotes the conditional expectation of X, given A. Let j be fixed. For 
each AT, set Zm — \Xj\ 2 if \Xj\ > e/\bN,j\, and otherwise. Since |&jv,j| — * 0, Zjy < 
\Xj\ 2 for each N, and Zm — ► almost surely, we may use the dominated convergence 
theorem to show that for each j = 1,2,. . ., EZ N = E(\Xj\ 2 | \Xj\ > s/\bisr,j\) — > as 
N -> oo. Thus, the third condition J]jLi ^(I^Vjl 2 | I^Vjl > e) ^ is satisfied. The 
conclusion now follows directly from application of the Lindeberg-Feller theorem. □ 

For each j = 1, . . . ,N — 1, we take bjq,i = and Xi = LUi in Proposi- 

tion and we see that Oj converges in distribution to a^ 00 " 1 as N — > oo, where 
a^ 00 ' is a Gaussian random variable. Moreover, due to the orthogonality of W^ N \ 
a\ N \ . . . , Im-Ii m some sense, become independent in the limit. 

Lemma 4.4. The random variables a\ , a 2 °° , . . . are independent. 

Proof. We need to show that for any finite A C N, the collection {a^JteA is a 
set of independent random variables. For simplicity, we only prove this for A — {1,2}, 
but a similar argument works for a general case. 

Let t\ and t 2 be given. Set 



VW+4 

Then, as N — > 00, max., \bpf j\ approaches zero because maxj \ W\j\ and max., \W2j\ go 

-3=1 °N,j 



to zero. Also, it is easy to check that Ylj=i b% ■ — 1 f° r au N. Applying Proposi- 



tion EH we see that 



+ „W -Lf n W N 

A 2 , / = E h »w -> ^ (°> 17 )> 

V t l+ t 2 j=l 
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which implies that tia[ N ^ +t2d^ — ► Af(0, o 2 (t\ + i§)); which in turn implies the con- 
vergence of the joint characteristic function of ^ and as N — > oo. Specifically, 

£:e itia5 oo) +it 2 a<° c ' ) = lim jBel (tia< JV) +t 2 a< N ») = e ~ CT 2 (t?+t| )/2 = g-o- 2 t?/2 g-^t=/2_ 



Therefore, aj 00 '' and ajj are independent. □ 



(oo) (oo) 



Let us come back to the proof of Theorem l4.2l As a consequence of a-j° ' , a£ , . 
being independent Gaussian random variables, we have 

q c < P({-eM < af } < -em, j = 1, . . . , N }) 



N 

-r 2 /9,T 2 , 

e 



em 
eM 



as N — > oo. Since this holds for any fixed No, and since the right-hand side goes to 
zero as No — > oo, we conclude that g c — > 0. This completes the proof of Theorem 14 .21 
□ 

5. Example. A network of voltage-controlled oscillator (VCO) devices can be 
built as an example of systems with diagonalizable interaction. The behavior of the 
jth VCO in the network is described by its phase variable 9j, which satisfies 

where ojj is the center frequency and Ij(t) is the input signal from other VCOs. The 
system is diagonalizable if, for example, Ij(t) has the form 

N-i / N 



= £ £<° sin 




k=i \e=i 

with WW defined by COl . 

This type of interaction can be implemented using commercially available circuit 
elements, as follows. The sine terms on the right-hand side can be constructed as 
sums and products of output voltages: 

\e=i / b t=i ^ 

where the sum is taken over all (ordered) binary AT-tuples b = (b\, . . . , 6jv), bg = 0, 1, 
such that Ylfbi 1S °dd- This means that I jit) is a sum of terms that are products 
of sin(W^9i) and cos(W^9f). Signals of the form sin(W^'9i) can be obtained 
by using the amplified version of the input Wcl Ii(t) as the controlling voltage in a 
separate VCO with center frequency W^'Ug. From these we can get coa(W^$i) 
by the phase shift of ir/2. Finally, I jit) is obtained by putting these signals through 
multipliers and adding the outputs. 

6. Discussion. In order to gain additional insights, let us consider the case when 
the functions pj = p do not depend on j, and have the range [—1, 1]. The arguments 
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Fig. 6.1. Probability of coherence g c (e; N) as a function of e for a = 1 and N = 100, 1000, 10000. 



in the proof of Theorem 14. 21 suggest that for a diagonalizable system (|l.l(l with large 
N, the probability q c that it is coherent is approximately 

N-l 



1c 



1 



V2 



7TO~ 



-i JV-1 









where erf(x) is the error function. Typical graphs of q c (e; N) are plotted in Fig. 16.11 
One can see that for any finite value of N, there seems to be a sharp transition point 
through which g c (e; N) changes from to 1. However, unlike the mean-field model of 
Kuramoto, this point keeps shifting to the right as N increases, and tends to oo in the 
limit of N — > oo, although it can be shown that this increase is at most O(yhiN): 
Lemma 6.1. Let a > and < q < 1 be fixed. Define e q ^(N) implicitly by 

<o; i.e., e n 



Then e q , a (N) = 0(VhiN) as N 



-q,cr 



(N)/VlnN 



q = q c (e q>a (N); N). 
bounded as N — > oo. 

Proof. Let x = (ay/2)- 1 e g><r (N). Then, using a known estimate for the error 
function, we have for x > 1 



erf (x) > 1 - 
Hence, we have the estimate 



2e- 



q > (1 - C Q e- x ) Ar - 1 > 1 - (JV - l)C e- x ~ 



> 1 - 



2e" 



Vtt + Vtt + 4 



where 



C 



/tt + V7T + 4 



Here we used the relation (1 — x) n > 1 — nx, which is valid for n > and < x < 1. 
The estimate for e q ^o-(N) can be obtained by rearranging: 

e q AN) < <7VCi+21n(iV-l), 
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where d = 21nC - ln(l - q). This implies e q . a {N) = 0(y/]nN). □ 

7. Conclusions. In this paper, we have defined a class of systems of phase oscil- 
lators characterized by having diagonalizable interactions. For a system in this class, 
complete separation of variables through appropriate changes of variable is possible, 
which enables us to draw rigorous conclusions about the probabilistic properties of 
the system. In particular, we have shown that partial coherence of the system almost 
surely implies coherence and, in the limit of large system size, the system is almost 
surely incoherent. A major implication of our result is that, unlike the mean-field 
model of Kuramoto, diagonalizable systems cannot exhibit a sharp transition from 
incoherence to coherence. This provides some insight into what is necessary to see 
such a transition in a system of phase oscillators. 
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